^''S:¥r 


Technic#; 


'/ '> ^ 'V v , ' 

\ 'V 


A 1 v. A '^'' 


2632 

,>..• 1 -* ( . , 

" . ;' i . . ■ ■•■ - “ • V •'■< r , ,-v -• 

' i,., ) ■ • •..» _-v> " ’• /. 

December 1986 






*- .■ :. >•• • . ■ ■ •, ■ • * . " f - 
I 2 ‘*t 


Raymond L. Barger 
and Allen K. Dominek 




1 Tr?* A 

, , ■ J ** ' 

X ^ -y A -W ' 1 
• ■ *■ ■ . , - *>- — . 


{NAS A -TP-2 6 32) TEEOEY FOE CC MEUTIN G THE 
FIELD SCATTERED FEOM A SNOCIH INFLECTED 
SOSFACii (NASA) 23 p CSCL 20F 


N87-13264 


fJrjclas 

44717 


NASA 

Technical 

Paper 

2632 

1986 


Theory for Computing the 
Field Scattered From a 
Smooth Inflected Surface 


Raymond L. Barger 

Langley Research Center 
Hampton, Virginia 

Allen K. Dominek 

Ohio State University 
Columbus, Ohio 


NASA 

National Aeronautics 
and Space Administration 


Scientific and Technical 
Information Branch 


Summary 

A theory is described for computing the reflected or scattered field from a smooth body 
with inflection points. These inflections occur in certain directions at each surface point for 
which the total (Gaussian) curvature is zero or negative. For surface illumination in one 
of these critical directions, the usual reflection formulas obtained by the high-frequency 
approximation are inapplicable, and a shadow zone exists in the reflected field. Scattering 
into the shadow zone is treated, as well as specular reflection. This theory should have a 
variety of applications such as for certain optics problems, computer graphics modelling of 
three-dimensional shapes, and the design and analysis of specialized microwave reflector 
antennas. 


Introduction 

When a wave is incident on a smooth surface, the reflected field can usually be computed, 
in the high-frequency approximation, provided that the curvature distribution of the surface 
in the vicinity of the reflection point is known. For the monostatic case (receiver coincident 
with source), the two principal curvatures are required. For the bistatic case, the curvature 
in the direction of the incidence plane is required as well as the curvature in the orthogonal 
direction. If one of these curvatures is zero, the usual formula is inapplicable. This 
problem is caused by the vanishing of the quadratic term in the local phase distribution 
function. The next higher order nonvanishing term must be included in the surface current 
integration to obtain the high-frequency approximation for the scattered field. In this 
analysis, it is assumed that this term is the cubic term. 

A special case of inflection point scattering has been treated in reference 1. The analysis 
of reference 1 is limited to reflection from bodies of revolution in a system configured in 
such a way that the body surface equation is expressed in a coordinate system whose origin 
is located at the ray source. 

Such a restricted analysis cannot be applied to more general configurations. For a body 
of revolution, all the required surface geometric parameters are readily obtained from 
the equation of the meridian line, but for more general surfaces this simplification is not 
available. Even for an axisymmetric surface, it would not be a trivial matter to adapt the 
analysis of reference 1 to the situation in which the body is illuminated from an off-axis 
source, especially if, for the bistatic case, the incidence reflection plane does not intersect 
the body axis. 

The present analysis is not restricted to axisymmetric shapes or to a special coordinate 
system. However, the analysis is limited to the far-field condition, primarily for clarity 
in presentation of the concepts. Modification of the analysis to account for near-field 
conditions should be fairly straightforward. 


Symbols 

A, B , a, b 

Af[] 

c 

E 

E 


parameters defined by equations (37) 
through (40) 

Airy function 

proportionality constant used in derivation 
of equation (7) 

electric field vector 

component of E 


E,F,G 


metric coefficients evaluated at r 0 


e,/,g 

e 

e Ts 


e 0 

e 0» 


e ±s 

H 

H 

I 

la 

b 

Jp 

J 

j 

k. 

m 

N 

P i t P s 

R 

R 

r 

r 0 

S 

s 

4 

U , V 


»7.£ 


curvature parameters evaluated at r 0 
unit base vector 

unit vector with direction of vector projec- 
tion of I s onto tangent plane 

unit vector e . x (see fig. 3) 

unit vector in scattering plane orthogonal to 
Is (see fig. 3) 

= TV x 

unit vector orthogonal to incidence plane 
unit vector orthogonal to scattering plane 
magnetic field vector 
component of H 


ray direction vector 
direction of I s — I{ 

vector projection of onto tangent plane 

projection of surface current field onto 
polarization plane of scattered ray 

component of Jp 

= v^T 

wave number 
magnitude 

local surface normal direction vector 

projection operators defined by equa- 
tions (50) and (51) 

rotation matrix 


distance from r 0 in direction I s 
surface vector 

point of ray incidence on surface 
surface area 

distance along arc on surface 

variable defined by equation (42) 

independent surface variables 

local distance of surface from tangent plane 
at r 0 

local surface variable, varying in e ^ or 
direction, respectively 
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0 


e s 

e* 

K 

X 

V 

$ 


angle of incidence 
scattering angle 

angle between I a and N, cos 0* = N • I a 

normal curvature of surface 

direction of surface arc at r 0 as indicated by 
dv/du on arc 

permeability 

phase 

angle between e-p s and e„ 


OJ 

Subscripts: 

i 


s 

e 


u,v,£,ri 


,-L 


angle between I a and 

angle between I a and 
circular frequency 


incidence 

scattering 

component in direction defined by intersec- 
tion of incident plane with plane normal to 
incident ray 

partial differentiation with respect to u, v, £, 
and r], respectively 

parallel to or orthogonal to incidence plane, 
respectively 


Reflection Shadows Created by Surface Inflection 

Figure 1 shows the shape of a body of revolution whose meridian line contains two 
inflection points. Arrows indicating the incident beam direction and the distribution of 
reflection directions demonstrate the existence of a reflection shadow region both for the 
monostatic case (fig. 1(a)) and for the bistatic case (fig. 1(b)). 

Figure 2 shows the distribution of specular directions for this body. It is seen that when 
the incident ray direction is near the horizontal, there is only one specular point. When 
it is near vertical, there are three specular points. These two regions are separated by the 
inflection points. 

These considerations apply, for purposes of illustration, to the axisymmetric shape 
illustrated in figure 1. The following analysis is not limited to axisymmetric shapes, 
however, but applies to arbitrary smooth shapes that are representable as a vector function 
of two surface variables, r = r(u, v) . A useful class of such surfaces is described in 
reference 2. The analysis represents, in some sense, a generalization of the analysis of 
reference 1. Another approach might be to treat the inflection condition as approximately 
representable as a near confluence of two reflection points and develop the analysis as in 
reference 3. 
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Analysis 


Local Surface Geometry 

Consider a body illuminated by a plane wave propagating in direction / 2 . The scattered 
field received by a far-held observer in direction I s is determined primarily by the local 
surface geometry near point r 0 at which the surface normal N satisfies the reflection 
condition 

Is = k - 2(? t - • N)N (1) 

If the surface is described as a vector function of two variables, r = r(u, v), then N is the 
direction of r u x r v . The local approximation to the surface at r 0 is 

1 2 

r = r 0 + (r u du + r w dv) + -(r U u dv 

2 1 o 

+ 2r uv du dv + r vv dv ) + — (r U uu du 

9 9 o 

+ 3r uuv du* dv + 3r U vv du dv * + r vvv dv 0 ) (2) 

The local shape of the surface is determined by its deviation from the tangent plane, that 
is, by the component of r — r 0 normal to the tangent plane, 


(r — r 0 ) • N = -( r uu ■ N du 2 + 2r uv ■ N du dv + r vv ■ N dv 2 ) + ■ • • (3a) 

The linear terms 

r u ■ N = r v ■ N — 0 

since r u and r v lie in the tangent plane. Equation (3a) can now be written as 


(r — r 0 ) ■ N = i(e du 2 + 2/ du dv + g dv 2 ) (3b) 

& 

where e, /, and g are the curvature parameters defined as 

e = r uu-N 
f — r uv • A 

g = *vv • A 

The incremental arc length is given by the formula (ref. 4, p. 58) 

ds 2 = E du 2 + 2 F dudv + G dv 2 (4) 

Here E, F, and G are the metric coefficients defined as 


E — Ty, ■ Yu 
F ~r u -r v 
G = Ty ■ Ty 

If a direction is specified on the surface at r 0 by fixing the slope A = dv/du, the normal 
curvature in direction A is 


(e + 2/A + g\ 2 ) du 2 e + 2/A + gX 2 
( ds/du) 2 du 2 “ E + 2FX + GX 2 


( 5 ) 
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The directions for which k ^ attains its extreme values are the roots of the quadratic 
equation 


A 2 

E 


-A 1 

F G 


= 0 


( 6 ) 


e f 9 


(See ref. 4, p. 80.) These are the principal curvature directions, and they are mutually 
orthogonal (ref. 4, p. 80). The maximum and minimum values of normal curvature obtained 
by substituting these roots into equation (5) are called the principal curvatures. 

The direction A,, in which the incidence plane intersects the surface at r 0 (see fig. 3) 

can be found by projecting the incident ray vector onto the tangent plane 



This vector is proportional to r u + r^A,, 


Ip — c(r u + r w A,,) 


Thus, dividing the scalar products 


J p -r u = c{ r 2 +r u ■ r^A,,) 


Ip ■ r v = c(r u ■ r v + r 2 AJ 


and solving for A,, yields 


A„ = 


(jp • 

r l ~ 1 

[ip • 

1 (r« • r v ) 

(ip ■ r u> 


{jp ■ ru j 

1 (r« • r v ) 


The direction A^ orthogonal to A |( is given by the formula (ref. 4, p. 59) 


( 7 ) 


( g + rA„) 

( F + GA ll) 

Now new local coordinates are defined by the relations 


or 


= dv — A^ du 


dr) — d£ 



(8) 

(9a) 

(9b) 
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(10b) 


A ± dr] - Ah d£ 



Denote the unit base vectors in the tangent plane in directions and A^, by ^ and 
e^, respectively. Thus, £ varies in the e„ direction and rj varies in the e ^ direction. To 


be more precise, £ varies along a curve formed by the intersection of the incidence plane 
with the surface, and r/ varies along a curve formed by the intersection of the plane of ft 
and with the surface. For these curves, the curvature vector at r 0 is in the ft direction 

when the curvature is not zero. If the curvature is zero, the vector d^r/ds^ points in the 
N direction. Arc length in the e,. direction is given by 




ds „ = 


E + 2FX U +G\l 


A„ - A 


di 


(Ha) 


Similarly, in the direction, 


Je + 2FX + GX 2 , 

ds ± = . _ » dr] (lib) 

A X A || 


Scattered Field Calculation 

According to reference 5 (p. 149), the field scattered in direction I s is obtained by 
integrating the physical optics surface current distribution: 


_ ~JU]l 

S 2t xR 


jk R J ft x H t - - ( 


N x H, ■ I, 




jkr(T s -T^j 


dS 


( 12 ) 


The phase function 

* = kr&r,) ■ (I, -li) (13) 

is expanded first in the e.. direction, taking the £, f] origin at the point of incidence, to give 


$ 

k 


r M + <■(( + + |j r £e £{ 3 


(h-?i) 


In the Cj| direction, the surface vector can be written, approximately, as 

r -r 0 (r?) = ds ||«|j 

where f is the distance of the surface from the tangent plane. Thus, 


(14) 


( 15 ) 
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f(o) = ^(o) = 0 



For simple monostatic or bistatic reflection 


I s — L = 2 cos 6 N 


Thus, for this case, the linear term is 

r £ ‘ (j s ~ A') d Z = 0 

The dot product in the quadratic term of the phase expansion becomes 

|r cr (? s -/ t -) = r tf -iVcos* 


( 16 ) 


(17) 


If r£c(£) is not identically zero, the remaining terms can be neglected, and the integration 
with respect to £ can be carried out by the stationary phase approximation. 

Now, considering the term r (rj) ■ ^ I s — in equation (14), similar considerations apply 

1 -. 2,1 


r(r]) — To + Tr) dr] + —rrjr] drj + ~r rjrfri dr) + . . . 


(18a) 

(18b) 


r{rj) = r 0 + ds e ± + $(r])N 
where f(0) = fy(0) = 0. Thus for simple reflection, 

*rr(ls-h) =0 (19) 

rtr) • (js ~ = r rm -N cos 6 (20) 

These expressions can be written in terms of the original surface variables u, v as follows: 


*u + A , r v 

*7) = fu u r) + Tv v r) — ~ r r 

A 1 “ II 


ds 


_L 


dr] 


— r T ] ' r T ] ~ 


E + 2FX l + GX\ 


A ± -A, 


2 2 Tuu 2 r m;Aj_ + r OT Aj^ 

Yrjri = Yuu u r] T 2ruv u r]Vr) + r vv^rj = 


A ± -A| 


( 21 ) 


( 22 ) 


(23) 


Now, by definition of the curvature coefficients, 


Yjp] ' N — 


e + 2fXj~ + gXj~ 

( A x - A ||) 2 


Tun ' N — tt - 


e + 2 /A . + gX 2 , 


tL ~j” It A ^ “h Cjr\ 


a i v u,/ / 


(24a) 


(24b) 
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Try?; ’ N — 


ds 


drj 


(24c) 


The quadratic term in the phase function (eq. (20)), therefore, becomes 

JL 


rrjr) ■ N cos 6r]" — k t cos 6 ( — 


ds 

drj 


V 


(25) 


Relations analogous to equations (21) through (25) could also be written for the variation 
in the e,, direction to yield 



or, locally near the stationary point, 

• (js - /,•) £ 2 = K|| cosOsj 


(26) 


Assuming that the dominant contribution to the integral of equation (12) comes from the 
stationary point r 0 and removing the surface current from the integral leave the following 
integral to be evaluated: 



where the stationary phase approximation has been used to evaluate both integrals. 


Now, the surface current is resolved into components: 

TV x Hj • g|| = (Nx H i± e ± ^ • = H i± (28a) 

N x Hj • = (^N x H^geg^j ■ e ^ = H^g sin (n, eg'j — —H^g cos 6 (28b) 

eg = e ± x ? t - (29) 

Equations (28) give the surface current as 

N x Hj = H{^e || — H t g cos 6 e ^ (30) 


For the reflection direction I s , the projection of N x H ? onto the polarization plane yields 
the following components: 


|| 


e 0 I eg = H i± cos 0 eg 


(31a) 
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(31b) 


(^Hjg cosde^ • = —H^cosde^ 

Substituting E field components in equations (31a) and (31b) by the relations 


H i± = —E ie 


UJ^l 

k 


Hj’fl l-' ; | 

w ujh 1 

and using the results together with equation (27) in equation (12) yield 


(32a) 

(32b) 


f M - e — ( E » ) 

\ Ej_ s ) 2 R || V ^i-L / 

This well-known result has been developed in this paper in this manner in order to 
demonstrate the analytical relationship of this result with that for inflection scattering 
which is given in the following section. 



Reflection at an Inflection Point 

The preceding analysis fails if either or vanishes, that is, if the surface is inflected 

at r 0 either in the direction of the incidence plane or orthogonal to it. According to 
equations (24a) and (24b) (or the equivalent equations for £ variation) this means that for 
A = A^ or A = A^, 

e T 2/A + gX ^ = 0 

This equation has solutions 

A _ -/ ± \/ / 2 ~ eg 
9 

Thus, if eg - / 2 > 0, equation (34) has no real root. The surface is convex or concave in 
all directions at r 0 , which is termed an elliptic point. If eg — / 2 = 0, equation (34) has one 
real root, A = —f /g, and r 0 is a parabolic point. The surface is inflected in one direction 
only. If eg — / 2 < 0, equation (35) has two real roots corresponding to two directions 
in which the surface is inflected. In this case, r 0 is a hyperbolic point. Directions A for 
which equation (35) is satisfied are called asymptotic directions. If a wave is incident at 
a hyperbolic point so that the incidence plane is aligned with either of the asymptotic 
directions, or is orthogonal to either of them, equation (33) fails. For a parabolic point, 
only two incidence plane directions present this problem. 

When the problem does occur, a shadow region is created in the reflected ray tube. 
This effect was mentioned in the first section and illustrated in figure 1. By differentiating 
the reflection condition (eq. (1)), 


(34) 

(35) 


-dls 

ds 



dN „ 

dF + 2 



N 


it is seen that an extreme value of I s is attained for 
dN/ds = 0, corresponding to the inflection condition. 


dl/ds = 0 which occurs when 
Thus, the ray reflected at the 
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(40b) 


2 T VV ' (js ~ = a 


where relation (24c) has been used in equation (40a). 

With this notation, the integral (corresponding to eq. (27)) to be evaluated is 


I j kr i T --tf dS = d8|| /“ * 

x r /("’ +te ±) * 

./ -00 


j AsS+Bs,, 

e V 11 11 / 


(41) 


The inner integral is reduced to the previous type by the simple change of variable: 


* b 

S± _ + 2 a 


(42) 


Equation (41) then yields the result 


Jtro- (/,-/,) c -tf /4a r°° t°° dal e >Hl +fls ii) 

J—oc J — oo 

= e J k r 0 (ls-h) e _ jb 2 /ia e _ jir/4 fn 2tt 

V « (3A) 1 / 3 


ds, 


L(3A)V3j 


(43a) 


Jkro- (Z-Z) e _ jb 2 /4a foe roo ^ 

J—oc J —oo 

= J kr o (h-Ii) jtf/Aa (-jir/A 


0S I ds* ± e 1 


,(M +Ss h) 


ds n 


/ 27r 1 

1 2?r \i 

B 

W kmt cos0*J 

1 (3A) X / 3 

_(3A)V 3 _ 


(43b) 


where the integrations have been evaluated by the usual formulas (ref. 6, p. 447). 
According to equations (28a) and (28b), the surface current is proportional to 


N x H, = (H i± e || - H ie cos 6 e ± j 


(44) 


Resolving the scattering direction vector I s into components in, and normal to, the tangent 
plane yields 

I s = sin 0 S zjs + cos @s N (45) 


where 


e Ts — cos 0s e,. + sin 4> s e , 


(46) 


Orthonormal field directions in the polarization plane associated with the scattering 
direction I s are 

e q s = — cos 0 S erp s + sin0 s N (47) 


tj_ s = —sm<ps e., + cos cps 
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Therefore, the component of the surface current field in the polarization plane is 

J P = N x Hi - (N x Hi • ? s ) 7 S 

= (~H t _ i_ sin 4> s — HiQ cos 9 cos (p s ) ej_ s 
— cos 9 S (cos <j> s Hi_ i_ — sin </> s cos 
k 

= — [ {-Erf sin <f> s + Eij_ cos # cos (p s ) e_L„ - cos 0 S (cos </> s E. q 

UJfi 

+ sin (f> s cos 6Eij _ ) eg s \ (49) 

where relations (32a) and (32b) have been utilized. The matrix operators P ? , P s , and R 
are defined as follows: 


P * s (i co°s«) ( 5 °) 

?) 01 ) 


/ COS 4>s 
\ - sin (ps 


sin (p s 
cos (ps 


(52) 


The operator P 2 projects from the incident polarization plane to the tangent plane. The 
operator R rotates in the tangent plane from the e,, direction to the ey s direction aligned 


with the scattering plane. Then the operator P s projects from the tangent plane to the 
polarization plane associated with the scattering direction I s . With equations (50) through 
(52), equation (49) becomes 


( JQs 

\J±s 


PsRP 


W/i 



(53) 


Substituting from equations (43b) and (53) into equation (12) yields the final result as 


( E 0s ^ — e Jtt/4 -ib 2 /Aa / k jkr 0 yl s 

\ E±s ) R V m/c _L cos 


x 1 \j B 

(3/A)i/3 [(3A) 1 / 3 

where R, A, b, and a are defined, respectively, by equations (37), (38), (39), and (40). If 
scattering only in the incidence-reflection plane is considered, 


PsRP, 


Eil) 


(54) 


sin 9* = 0 


consequently, 


6 = 0 

and the factor e ^ !^ a drops out of equation (54). In this case, the well-known inverse 
square-root dependence on the normal curvature k ^ is obtained. However, if scattering 
outside the incidence reflection plane is considered, 


e -jb 2 /Aa ^ 1 
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inflection point r 0 represents the edge of the shadow region. As was mentioned earlier, 
the field scattered into this shadow region is significant. In order to calculate it from the 
surface current integral (eq. (12)), the cubic term must be included in the expansion of 
the phase function in the inflection direction, and scattering directions different from the 
reflection direction must be considered. 

Assume that the surface is inflected at r 0 in the direction. Then, in the phase 

function expansion (eq. (14)), the quadratic term vanishes. Since the reflection condition 
is no longer assumed, equation (1) is no longer applicable. Let I a denote the direction of 
I s — I{ and m denote its magnitude. Then I s — I{ can be resolved into components: 

I s — /j = ml a = m (cos 0* N + cos^e^ + cos^y (36) 

Since, from equation (15), r^(o) = ey, the linear term in the phase function becomes 


• (l s — I^j d£ = km cos ^|| dsy 

(37a) 

kr ^ • (l s — I^j d£ = B day 

(37b) 

or locally, 


kr^ ■ (l s - l^j d£ = Bs|| 

(37c) 

Differentiating equation (15) three times and using equation (36) yields for the cubic 
term, 

kT iU ' (j s ~ A’) £ 3 = k ‘ mcos(? * £ 3 

(38a) 

kr ggg ■ (is - IjJ £ 3 = k-~-m cos 0* s| 

(38b) 

kr tt£ ■ (js ~ %) £ 3 = As| 

(38c) 

Now, considering the phase variation in the direction, equation (19) 

is replaced by 

**■(/„ l)=kmcos4, ± 

(39a) 


(39b) 

Equation (20) for the quadratic term is replaced by 


-r vr] ■ (l s ~l{j = —$r)r]{0)N • (cos 9* N^j 


= cos 8* ■ N 

2 ds\ 

km .* 

= — — — K cos 9* 

(40a) 
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and an additional dependence on k ^ occurs through the parameter a. (See eqs. (40a) and 
(40b).) 

Furthermore, for the reflected ray that satisfies the reflection law (eq. (1)), 

% = N 


consequently, 



= 0 


in equation (36) and therefore, B = 0 (eqs. (37a) and (37b)). In this case, the factor 
involving the Airy function in equation (54) becomes 


1 

(3A) 1 / 3 


A*[0] = 


1 

3A 1 / 3 T(2/3) 


(55) 


Now consider the case for which the surface is inflected in the direction. An 
analysis similar to the preceding one leads to a result like equation (54) with £, r) variables 
interchanged and replaced by /c... A difference in the physical phenomena, however, 


should be noted. For this case, the reflection shadow occurs because of the inflection in 
the e ^ direction, and consequently the incidence- reflection plane is the boundary of the 
shadow region. In order to compute scattering into this shadow zone, one must treat those 
scattering directions for which b ^ 0, that is, for which the factor /Aa occurs j n the 
field equation. For the previous case, on the other hand, one can study the diffraction into 
the shadow region for scattering directions in the incidence-reflection plane. 


Location of Inflection Points 

For the far-field monostatic problems, specular points are normally located by compar- 
ing the local surface normal vector N with the reflection direction 7 r . For the corresponding 
bistatic problem, N is compared with the direction of I r — I t ■ To determine whether a specu- 
lar point is elliptic, parabolic, or hyperbolic, the sign of the total curvature (or equivalently, 
the quantity eg — / 2 ) is required. For the general scattering problem, it is advantageous to 
know the distribution of parabolic and hyperbolic points on the surface. This distribution 
can be determined by making an orderly survey of the quantity eg = — / 2 over the surface. 
Its zeros determine the lines of parabolic points which separate the hyperbolic regions from 
the elliptic regions. 

Figure 4 shows a surface having both elliptic and hyperbolic regions, with the line of 
parabolic points separating them. 

Concluding Remarks 

A theory has been presented for computing the reflected or scattered field from a smooth 
body with inflection points. Although the analysis was developed for electromagnetic 
waves, the results can easily be reduced to apply to the acoustic problem. Far-field 
conditions were assumed. The only restriction regarding surface geometry, other than 
its smoothness, was that it can be represented, at least locally, as a vector function of two 
variables. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
August 15, 1986 
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Shadow boundary ray 
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(a) Monostatic. 

Figure 1. Diagram indicating reflection shadow formation at inflection point. 




Figure 1. Concluded. 
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Figure 3. Geometry at point of ray incidence r 0 . 
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